Oxygen Vacancy Induced Ferromagnetism in Y 2 0^- x 
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Ab initio calculations within density functional theory with generalized gradient approximation 
have been performed to study the effects of oxygen vacancies on the electronic structure and mag- 
netism in undoped V^Os-rc (0 < x < 0.5). It is found that the introduction of oxygen vacancies 
would induce ferromagnetism in V205_a; with the magnetization being proportional to the O va- 
cancy concentration x. The calculated electronic structure reveals that the valence electrons released 
by the introduction of oxygen vacancies would occupy mainly the neighboring V d^y-dominant band 
which then becomes spin-polarized due to intra-atomic exchange interaction, thereby giving rise to 
the half-metallic ferromagnetism. 
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Ferromagnetism in diluted magnetic semiconduc- 
tors (DMSs) has recently received a lot of attention 
for fundamental research and possible applications in 
spintronics [j], 0, H|. Although a number of transition- 
metal (TM) oxide based DMSs, such as Co:Ti02 and 
Co:ZnO, have been reported, 0, @ the origin of ferromag- 
netism in many cases is still under much debate. [fj, |7|, [8| 
The controversy of the magnetic origin lies in the struc- 
ture of DMSs in the first place, i.e., whether it is due to 
the intrinsic contribution from the DMS phase or from 
some secondary phases. [IJ It is noticed that even if the 
formation of the secondary phases could be ruled out, the 
valence effect on ferromagnetism, the role of the micro- 
scopic complex structure formed by oxygen, dopant, and 
host TM elements need to be clarified. Oxygen vacancies 
may also play an important role for the ferromagnetism 
in oxide based DMSs because oxygen vacancies not only 
modulate the valence of neighboring TM elements but 
also cause a change of the band structure of host oxides; 
both factors in turn can make significant contribution to 
the ferromagnetism. Q Several conflicting experimental 
results of oxide based DMSs could be partly related to 
the difference in their oxygen vacancies because of dis- 
tinct sample growth conditions. 0, Q . Therefore, under- 
standing of the role of oxygen vacancies in TM oxides 
is essential to clarify the origin of the ferromagnetism in 
the oxide based DMSs. 

Vanadium oxides consist of compounds with single 
[V 2 5 (V+ 5 ), V 2 4 (V+ 4 ), V 2 3 (V+ 3 ) and VO (V+ 2 )] 
and mixed valences which form a fascinating class of 
multifunctional materials. Vanadium oxides have been 
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widely used in catalytic, optical, electrical and electro- 
chemical devices. [13, EH, EH V2O5 is a diamagnetic in- 
sulator at room temperature with an orthorhombic lay- 
ered structure. It undergoes a phase transformation to a 
metallic phase at about 530 K.pjl, [l4| By comparison, 
V2O4 is a diamagnetic semiconductor with a tetrago- 
nal rutile structure below a critical temperature of ~340 
K and a metal with a monoclinic structure above. [l5| . 
While the physical properties for various vanadium ox- 
ides have been extensively studied, [l(| E3, EH the role of 
oxygen vacancies on the structural, electronic and mag- 
netic properties of V2O5 remains yet to be explored. We 
have therefore performed ab initio theoretical studies of 
the electronic structure and magnetism in V2G , 5_ ;r . Our 
calculations reveal that the presence of oxygen vacancies 
introduces electrons onto the neighboring V d xy domi- 
nant conduction band which is then spin splitted due to 
intraatomic exchange interaction, leading to the forma- 
tion of the magnetic order. Our finding calls for critical 
re-examinations of the role of the intrinsic O vacancies in 
the formation of magnetism in TM based DMSs, and also 
suggests a new direction for search for ferromagnetism in 
semiconductors . 

Our ab initio calculations are based on density func- 
tional theory with the generalized gradient approxima- 
tion to the electron exchange-correlation interaction. . 
Two approaches are used to model the oxygen vacan- 
cies in V205„2;. One is the so-called virtual crystal ap- 
proximation (VCA). In V205_a;, x O atoms are removed, 
and this is equivalent to increase the number of the con- 
duction electrons by 2x. In the VCA, this is simulated 
by artificially increasing the electron and atomic num- 
ber of V from 23.0 to 23.0 + x, i.e., replacing the V 
atoms by the virtual atoms with an atomic number of 
23.0 + x. We have performed self-consistent band struc- 



FIG. 1: (Color) Isosurface of the spin density of 
(x = 0.125) in the 1x2x2 supercell. The numbers are the 
numberings of the V atoms in the unit cell that are refered to 
in the text and Table I. 



ture calculations for V2O5-X with several x values in 
both nonmagnetic (NM) and ferromagnetic (FM) states. 
We use the highly accurate full-potential linearized aug- 
mented plane wave (FLAPW) method, as implemented 
in WIEN2k package. [2(| The experimental lattice con- 
stants used here are a — 11.512 A, b = 3.564 A, and 
c = 4.368 A. [2l| The cutoff angular momentum for the 
spherical harmonics expansion is 10 and as many as ^120 
augmented plane waves per atom are included. 

The other approach is the supercell (SC) method which 
is much more CPU time-consuming especially for small x 
values. Nonetheless, the SC method is more realistic. In 
particular, we can study local structural, electronic and 
magnetic properties in the vicinity of the O vacancies. 
We have carried the SC calculations for the x = 0.25 and 
0.125 cases using the 1x1x2 and 1x2x2 supercells, 
respectively. Pure V2O5 has a layered structure with two 
formula units (f.u.) per unit cell. [2lT| Each layer contains 
edge- and corner-shared VO5 pyramids with three types 
of O atoms, namely, apical oxygen (01), chain oxygen 
(02) and bridge oxygen (03), as illustrated in Fig. 1. 
In order to investigate which kind of the O vacancies is 
energetically favorable, we perform ab initio structural 
optimizations using the 1x1x2 supercell with one 01 
or 02 or 03 vacancy, respectively. We use the computa- 
tionally more efficient projector augmented wave (PAW) 
method. [22l| as implemented in VASP package. [23| A 
large cutoff energy of 400 eV is used, and the equilibrium 
atomic positions are obtained when the forces acting on 
all the atoms were less than 0.02 eV/A. Our 1 x 1 x 2 SC 
calculations show that the total energy per vacancy for 
the 01 vacancy is 1.41 and 2.15 eV lower than that of the 
02 and 03 vacancies, respectively, suggesting that apical 
oxygen (01) is the preferred O vacancy site. Thus, we 
perform the structural optimization for the 1x2x2 SC 
with only the 01 vacancy. Our 1x1x2 and 1 x 2 x 2 SC 
calculations show that the 01 vacancy in the FM state 
has a lower total energy than in the NM state, respsec- 
tively, by 78 and 86 meV per vacancy, showing that the 
ferromagnetism is stable in V2 05_a;. The stable FM state 
is also found in the VCA calculations for all the x values 
considered. We also perform the self-consistent spin po- 
larized calculations for the antiferromagnetic (AF) con- 
figuration with the initial magnetic moment on the V 
atom next to the vacancy being antiparallel to the other 
V atoms. However, the self-consistent calculations al- 
ways converge to the FM state, indicating that the AF 
state is unstable. 

Final self-consistent band structures for the 1x2x2 
supercell with the theoretically determined atomic struc- 
ture were calculated using the FLAPW method. The 
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FIG. 2: Total and site-decomposed density of states for 
V2O5-X with (a) x = 0.0, (b) x = 0.125 (SC calculation) 
and (c) x = 0.25 (VCA calculation). 



calculated V local magnetic moments for x — 0.125 are 
listed in Table I. The total magnetic moment is 2.0 [is 
per vacancy, and is the same as that obtained from the 
VCA calculations. Majority of the magnetization is lo- 
cated on the V atoms (d^-orbital) next to the O vacancy 
(V7,V3,V9 in Fig.l) and also in the surrounding intersti- 
tial region, as can also be seen from Fig. 1 which shows 
the spin density distribution. The induced magnetiza- 
tion on all the O atoms are negative (i.e., antiparallel 
to that on the V atoms) and the sum of the O atomic 
moments is -0.360 fiB per supercell. However, the lo- 
cal magnetic moments on most O atoms are negligible 
(< 0.01/is), though six O atoms have a magnetic mo- 
ment of 0.01 - 0.05^ B . 

Fig. 2(a) shows that pure V2O5 is a semiconductor 
with a band gap of ~1.8 eV. While the valence band is 
clearly O p dominant, the lowest conduction band is of 
mainly V d xy character. Interestingly, Fig. 2(b) shows 
that when one O vacancy is introduced, the lowest con- 
duction band is strongly spin-polarized, and the spin up 
one is partially occupied whilst the spin down one is still 
empty, resulting in a FM state with a spin moment of 2 
fiB per vacancy. Clearly, when one O atom is removed, 
the two electrons released are added to the narrow V 
d xy dominant conduction band. However, because of the 
rather strong V intra-atomic exchange interaction (I ex ) 
of ^0.354 eV,[24[ this partially occupied V d band be- 
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TABLE I: Local magnetic moments (m s ) (hb /atom) of V 
atoms in the ferromagnetic V2O5-X (x = 0.125) from the 
1x2x2 supercell calculation. The total magnetic moment 
is 2.0 hb per O vacancy. The locations of various V atoms 
are indicated in Fig. 1. Int. denotes the magnetization in 
the interstitial region. Also listed are the V-site decomposed 
density of states at the Fermi level N(Ef) (states/eV/atom) 
in the nonmagnetic V^Os-a;. 
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FIG. 3: Magnetic moment per V atom as a function of O 
vacancy concentration x from VCA and SC calculations. 



comes magnetically unstable and spin split, giving rise to 
a ferromagnetic order, as the spin-polarized band struc- 
ture for x = 0.125 and 0.25 in Fig. 2(b-c) shows. To 



corroborate this Stoner mechanism for ferromagnetism, 
we also calculate total and site-decomposed density of 
states in the nonmagnetic A^Os-^ (x — 0.125) in the 
1x2x2 supercell. Table I shows that V7 next to the 
O vacancy (see Fig. 1) has such a large local density of 
states at the Fermi level (N(Ep)) that the Stoner cri- 
terion (I ex N{Ep) > 1) for the magnetization formation 
is satisfied. Fig. 3 shows the calculated magnetic mo- 
ment per V atom as function of O vacancy concentra- 
tion (x). Remarkably, the calculated magnetic moment 
per V atom are proportional to x. This is due to the 
fact that the additional 2x electrons released by remov- 
ing x O vacancies go exclusively to occupy the lowest 
V ci-conduction band of the majority spin, as mentioned 
above. 

Summarizing, our ab initio calculations show that 
V2O5-X, with 0.0 < x < 0.5 is ferromagnetic. The two 
valence electrons released by the introduction of each 
oxygen vacancy would occupy mostly the d xy orbital 
of the V atoms next to the vacancy, and become fully 
spin-polarized due to intra-atomic exchange interaction, 
giving rise to the ferromagnetic order with a magnetic 
moment of 2/is per vacancy. We note that recent ex- 
periments indicate possible room-temperature ferromag- 
netism in oxygen-deficient V2 05-Z [25|]. We also note 
that the recent ab initio calculations indicate the un- 
expected ferromagnetism observed in HfC>2 being also 
driven by intrinsic O vacancies. [2g Nevertheless, the mi- 
croscopic mechanisms in the present and previous cases 
would be fundamentally different since in Ref. (26|, the 
magnetic moment is centered at the O vacancy. Clearly, 
the present work calls for critical re-examinations of the 
magnetism in TM oxide based DMSs especially the role 
of the intrinsic O vacancies. The present work also sug- 
gests a new direction for search for ferromagnetism in 
semiconductors . 
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